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Abstract. 

In his comment, Riemann defends the conventional kinetic Bohm criterion on the 
basis that the underlying approximations become rigorous in the limit Xd/1 0, 
where Ad is the Deybe length and I is a collision length scale for the dominant 
collision process. Here, wc expand on our previous arguments showing that the basic 
assumptions for the justified use of the conventional kinetic Bohm criterion are typically 
not satisfied in laboratory plasmas. We corroborate our argument with experimental 
data, as well as data from numerical simulations, showing that the conventional 
criterion is violated in common situations. In contrast, a formulation based on positive 
velocity moments of the kinetic equation provides a criterion that both agrees with the 
experimental data and reduces to the traditional expectations from fluid theory in the 
appropriate limit. 
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In his comment [T], Riemann defends the conventional kinetic Bohm criterion 



using theoretical arguments based on the Vlasov equation and the assumption of a 
collisionless, or infinitely thin, sheath Xd/1 0. Here is the Debye length, which 
is the characteristic scale of the sheath, and / is the characteristic scale of the plasma 
(or presheath), which is typically the shortest ion collision mean free path. Such two- 
scale analyses have been very successful at modeling important quantities such as ion 
and electron density, electrostatic potential and current profiles in the plasma-boundary 
region when Ad// ^ 1 (SJIUIHIIHIITIIHIIHI HO]. Some of these also provided models of 
the ion velocity distribution function (IVDF) [3l El El |9l [10] that have been shown to 
quahtatively predict experimental [HI [T2] and numerical simulation data [131 [IH HSl [IS] • 
However, there is one important difference; the IVDF is always empty for Vz > at the 
sheath edge in the theories, but not in the data. Theoretically, this is a consequence of 
the infinitely thin sheath assumption (no ions can escape the sheath). Experimentally, 
the sheath has a finite thickness, and a small fraction can escape either from scattering 
or from ions born from ionization that have sufficient directed energy to escape the 
sheath. Although this fraction is small, it has major consequences for the applicability 
of Eq. ([T]) because the left side diverges in this situation. In this Reply, we reiterate the 
main point of Ref. [IT]; Eq. ([T]) is not valid for arbitrary distribution functions, including 
cases experimentally measured. Instead, an alternative KBC based on positive-exponent 
velocity moments should be used. 

Bohm's original criterion [l8], Vi > Cs, was derived using two-fluid theory in 
the cold ion limit. The goal of a kinetic Bohm criterion is to generalize the fluid 
approach to account for arbitrary electron and ion distribution functions. Equation ([T]) 
applies to only a subset of possible ion distributions because the derivation assumes 
fi{vz = 0) = 0. Although this condition is satisfied for theoretical IVDFs calculated 
in the limit A/?// — ?■ 0, it is not typically satisfied experimentally because sheaths have 
a finite thickness. Equation ([T]) also assumes the distribution functions are Vlasov 
solutions, and utilizes properties such as fi is a function of only. However, non- 
Vlasov distribution functions are interesting because the presheath electric field always 
generates a current, implying that ions and electrons are not both Vlasov solutions in 
the same reference frame. One must be cognizant of the errors introduced when using 
non- Vlasov solutions in a theory based on the Vlasov equation. 

We are not aware of an experiment designed specifically to test Eq. ([T|, but 
experimental and numerical data has appeared in the literature that appear to violate 
it. Several experimental papers have presented measurements of the IVDF through 
the presheath and sheath regions, including the sheath edge. The primary experimental 
technique has been laser induced fluorescence (LIE) [IIl|12l[l9l|2ni[2lll22l[23l[2ll|25l|26l 
[271 1211 [29]. This has been applied in a variety of situations including single ion species 
plasmas [IIl[I21[ISl[2ni[2Il[221[22l[2a[251l2S], muhiple ion species plasmas [251 12Z], 
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and RF plasmas [28l 129] . A variety of numerical simulation techniques have also been 
applied to calculate the IVDF in the plasma-boundary region |13l [15], [161 iSOl [31] . IVDFs 
have also been measured near double layers, which are sheath-like structures where the 
Bohm criterion also applies, using LIF [32] and retarding field energy analyzers (RFEA) 
[55] , as well as calculation from particle-in-cell (PIC) simulations [311I3S]- Each of these 
references provide example IVDFs with a finite number of particles for t;^ > at the 
sheath edge. 

Three of these examples are shown in Fig. [T} Panel (a) shows experimental LIF 
data of an argon plasma from [26]. Panel (b) shows data from a hybrid simulation using 
a PIC routine for ions, and assuming electrons obey the Boltzmann density relation [12] . 
Panel (c) shows data from a PIC simulation pushing both ions and electrons [3T]. In 
each of these examples, there is a nonzero contribution to the IVDF for f ^ = at the 
sheath edgejf] Consequently, the left side of Eq. ([T]) diverges and the resultant criteria 
is not valid. It is useful to note that although the Xd/I limit is used to identify the 
sheath edge in the conventional theory, the sheath edge is experimentally meaningful 
even when this quantity is finite, so long as it is small. A variety of techniques are 
availible for estimating the location of the sheath edge. In fact, for the three sets 
of data in Fig. [l| three different techniques have been used. The experiments in (a) 
used a technique based on emissive probe theory [SH], estimating the sheath edge at 
z = 6.0 ± 0.5 mm. The hybrid simulations in (b) used the location for divergence of 
E from the cold Tonks-Langmuir problem, estimating the sheath edge at x/L = 0.917. 
The full PIC simulations in (c) used an estimate based on a coUisionless approximation, 
taking the sheath edge to be where the electrostatic potential has dropped Te/2 from 
the bulk plasma value; leading to an estimate of dX^ from the boundary. In each case 
a meaningful sheath edge can be identified, but Eq. ([T]) does not provide a meaningful 
Bohm criterion since the left side diverges. 

Reference [17] considers a couple of possible reasons for this disagreement. One 
of these is related to a step in the derivation of Eq. ([T| where integration-by-parts is 
used (see section 2.3 of [T7] and section 3c of [37]). Since this step is only valid for 
continuously differentiable functions, the restriction fi{vz = 0) = must be enforced. 
Instead of applying the integration-by-parts, and thereby restricting the set of IVDFs 
the theory can be applied to, one may write the criterion in the more primitive form 
from before this step is taken 



Equation ([2]) should apply to an expanded set of distribution functions compared to 
Eq. ([T|. However, the derivation of Eq. ^ still relies on symmetry properties of Vlasov 
solutions, and the errors introduced by using non- Vlasov solutions in this equation may 

I It is hard to tell what the ion density is at = for the LIF data in case (a) because the signal 
in this region of velocity space is within the noise. However, it is also hard to apply Eq. ([T]) to this 
data because the v~'^ moment amplifies the noise. For an example measurement with a clearer signal 
at i;^ = see Fig. 3 of [24]. 




(2) 




Figure 1. IVDFs in the plasma-boundary region from (a) LIF measurements from 
[26], (b) Hybrid simulations from [15 , and (c) PIC simulations from [31]. Reprinted 
with permission from: (a) G. D. Severn, X. Wang, E. Ko and N. Hershkowitz, Phys. 
Rev. Lett. 90, 145001 (2003). Copyright 2003, American Physical Society, (b) T. 
E. Sheridan, Phys. Plasmas 8, 4240 (2001). Copyright 2001, American Institute of 
Physics, (c) W. J. Miloch, N. Gulbrandsen, L. N. Mishra and A. Fredriksen, Phys. 
Plasmas 18, 083502 (2011). Copyright 2011, American Institute of Physics. 



be significant. Although Eq. ^ arises as an intermediate step in the derivation of 
Eq. ([T]) , Riemann argues [U EHl EH] that there is no useful sheath-related information in 
Eq. ([2]) because it is satisfied everywhere and that the integration-by-parts step, with the 
corresponding fi{vz = 0) = restriction, is required to obtain a meaningful criterion. 

For a specific example, consider the predictions of Eqs. ([!]) and ^ in the 
two- fluid limit, where electrons and ions have Maxwellian distributions: fe = 
neexp(— t;^/f|,g)/(7r^/^f|.g) and fi = ?T,jexp[— (v — Vj)^/ff,-]/(7r^/^f|,-) where Vj = ViZ. 
For these, Eqs. ([T]) and ^ reduce to 



1 . e"^' T,- 



h{V,/vT^) = ^ / dx- — < (3) 



00 



h{V,lvT^) = / dx "^f/, < ^ (4) 



and 

J~oo X-Vi/vTi Te 

respectively. We evaluate these integrals numerically by considering the Cauchy 
principal value: = \im^^o{f^'^^' ^ + f^^^^__^J- Figure 2 shows that as e gets 
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Figure 2. Plots of Ii and I2 from Eqs. (|3| and Q evaluated numerically excluding 
the interval Vi/vni^- The integral I2 converges as e — 0, but the integral Ii diverges. 
Two values of the ion flow speed are shown for each integral: Vi/vTi = 1 and 2. 




Vi/vTi 

Figure 3. The integral I2 as a function of Vi/vTi (solid line). Also shown are the small 
argument (dashed line) and large argument (dash-dotted line) asymptotic expansions. 



small, the integral I2 converges, but Ii diverges. So, in the two-fluid limit, Eq. ([T]) 
reduces to 00 < Ti/Te. Riemann argues in section 4.1 of Ref. gj that Eq. ^ can return 
fluid results by expanding [vz — V^)"^ for 3> f ^ in the integrand of Eq. ([s]). Indeed 
this does return Bohm's Vi> Cs criterion. However, the result is fortuitous because the 
full integral actually diverges independent of the magnitude of Vi. Such an expansion 
neglects the low velocity ions argued to introduce important new physics in the kinetic 
formulation [.Ij, but which cause the integral to diverge. 

Unlike Ji, the Cauchy principal value of the integral I2 converges and has an analytic 
solution 

/2 = -1 + ^{V,/vT^)e-''"/^T.evf,{V,/vT^), (5) 
where erfi is the imaginary error function (this is real for a real argument Vi/vTi)- 
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A plot of I2 is shown in Fig. ([s]), along with the large argument expansion [I2 = 
v\J{2V^) + ?)v\.J{W^)\ and small argument expansion il-i = — 1 + 2Vj^/f|,-). For 
^j/^e ^ 0.285, there are two solutions of Eq. Q. In the large argument limit, it 
reduces to Bohm's criterion > c^. In the small argument limit, it reduces to 
< VTiiX + Ti/Tf,) / \/2. The latter prediction is at odds with the fluid theory 
formulation, as well as experimental results. It is also noteworthy that for Tj/Te > 0.285, 
Eq. ^ is satisfied for any ion flow speed; there is little sheath-related information 
in this equation, as Riemann has shown previously from basic theoretical arguments 
[381 El]. Thus, in the fluid limit, we find that the left side of Eq. ([T]) diverges leading 
to an incorrect prediction. This divergence can be avoided by using the more primitive 
Eq. ([2]), but there is little sheath- related information in this equation. 

Reference [T7] considers another possible source of disagreement between 
experiments and Eq. ([T]), related to errors associated with using non-Vlasov distributions 
m Eq. Q. One way to justify the theory would be to include the collision terms in the 
analysis, then demonstrate that they are small in the final expression]!] As discussed in 
[1] and [17], keeping collisions (5) leads to the additional term 

on the right side of Eq. (JT]). In his comment, Riemann points out that if S is an even 
function of Vz-, this term will vanish because the positive and negative velocity intervals 
will exactly cancel. In this way, the Vlasov formulation may be justified for symmetric 
collision terms. The seminal theories of Tonks and Langmuir [S] , Emmert et al. [8J and 
Bissell and Johnson [9] are all two-scale theories that apply symmetric collision terms. 
These all lead to model distribution functions that satisfy Eq. (JT]). Riemann also points 
out that in the limit \d/1 0, theories that obey Eq. (JT]) can also be constructed 
using asymmetric collision terms because — t- 00 at the sheath edge in this theoretical 
limit [10]. 

Although theoretical models satisfying Eq. (JT]) can be constructed based on 
symmetric collision terms, or asymmetric collision terms in the \d/1 limit, weakly 
collisional plasmas in the laboratory have both asymmetric collision terms and finite 
Xd/I- The parity of ionization source terms is typically determined by the background 
neutral gas, which is often Maxwellian and leads to a symmetric source. However, sink 
terms are typically asymmetric in the laboratory frame. A frequently occurring example 
is ion charge exchange collisions where ions are removed from a flowing distribution. 
Another example is the sink associated with electron-neutral impact ionization. Near 
the sheath edge, the tail of the EVDF is often depleted in the direction facing away 
from the boundary because this population is lost to the boundary. Since tail electrons 
are responsible for the majority of ionization, this leads to an asymmetric sink term. 

§ Here "collision term" refers to all terms on the right side of the kinetic equation, including ionization 
sources and sinks, collisions of charged particles with neutrals. Coulomb collisions, and wave-particle 
scattering. 
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Ion-electron Coulomb collisions also generate an asymmetric collision term because 
the ion flow is shifted in reference to the electron distribution. This is a ubiquitous 
example because the presheath electric fleld always drives flow, leading to a nonzero, 
and asymmetric, collision term. Wave-particle scattering from flow-driven instabilities 
is another example that generates an asymmetric collision term |3T]. 

Because Eq. ([T]) leads to divergent integrals and places undue importance on 
ions with small velocity in the z direction, we suggested in Ref. [T7] that the theory 
should be developed from positive exponent velocity moments instead. There is also 
a semantic reason for using positive exponent moments. Bohm's original criterion [18] 
is a condition pertaining to the ion fluid flow speed at the sheath edge: Vi > Cg. A 
kinetic generalization should pertain to a kinetic analog of the ion fluid flow speed. The 
f ~^ moment of fi does not have a clear connection to the fluid flow speed. It is well 
known that a hierarchy of fluid equations can be rigorously developed from positive 
exponent velocity moments of the kinetic equation. With this approach, the fluid flow 
velocity naturally arises through the moment expression = J^^d^v\fs/ns, where 
the species s density is related to fs by Ug = J^^d^v fs- Reference [17] works from 
the philosophy that a kinetic generalization of Bohm's criterion should pertain to this 
kinetic analog of the fluid flow velocity. This approach provides an obvious physical 
correspondence between the kinetic and fluid theories, and simply returns the two-fluid 
results in the limit of thermodynamic equilibrium amongst individual species; the limit 
for which Eq. ([T]) was just shown to diverge. 

Reference ^7] shows that using the flrst two fluid moments of the full kinetic 
equation, along with Riemann's sheath criterion provides the relation 



Here density and flow velocity are deflned in terms of the previously mentioned 
deflnitions, and the other deflnitions are the temperature Tg = J^^d^vmsV^fs/{3ns) = 
stress tensor II^ = J^^d^v nigiyr^r — w^X/3)/s, and friction force density 
Rs = f^^d^v nis^^S^fs). Here v,. = v — and S{fs) can be any source, sink or 
collision term. Although Eq. ([T]) has an obvious connection to fluid theory, it also has a 
kinetic interpretation and can be written explicitly in terms of the distribution functions 
by substituting the moment deflnitions for the fluid variables. One difficulty with this 
approach is that it depends on spatial derivatives of temperature and stress moments. 
Most of the time, but not always (e.g., see section 4.2.2 of [E]), the terms in parenthesis 
in Eq. ([T]) are small because E becomes large at the sheath edge. This is related to the 
A^j/Z ^ 1 ordering. In difficult cases, higher order moments and a closure scheme can 
be applied to deal with gradient terms. 

Finally, we use the numerical data from Fig. [T|d to test the two criteria from Eqs. ([T]) 
and ([T]). Considering the = O.OlTe data from Sheridan's simulations. Fig. |4] shows 
the density, flow speed and temperature calculated from the appropriate moments of 
the IVDFs shown in Fig. [T|3. For this situation (with Boltzmann electrons), Eq. ([T]) 




qsUs - {risdTs/dz + dllzz,s/dz - Rz,s)/E 
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Figure 4. Moment parameters calculated from the IVDFs shown in Fig. [Tjo for the 
Tg = O.OlTe case. The sheath edge was predicted by Sheridan [T^ to be 0.917. The 
star shows the predicted speed at the sheath edge using Eq. Q. 



reduces to > \/ + v^i/"^, leading to the prediction Vi = Cs\/l + Tf,/Ti ~ 1.04cs at 
the sheath edge. This prediction is represented by the star in the figure where the sheath 
edge, as determined from cold Tonks-Langmuir theory by Sheridan, is predicted to be 
x/L = 0.917. As Fig. |4] shows, this prediction agrees well with the numerical results. 
Conversely, if we try to apply this data to Eq. ([T]), we find that J d^v fi/v1 — ?■ oo. 
Taking /g to be Maxwellian, which is consistent with the Boltzmann relation used in 
the simulations, the right side is ne/Tg. Thus, Eq. ([T]) provides the criterion oo < rie/Te, 
which is not consistent with the data. Applying the same analysis to the other data in 
Fig. [l] similarly leads to divergences in Eq. ([T]) , but the fiuid moment approach predicts 
Vi > Cs, in agreement with the data. 

In conclusion, the conventional KBC places undue importance on low velocity ions, 
which leads to divergences when the ion distribution function has particles at Vz = 0. 
We find that it is more productive to work from the perspective of generalizing the fiuid 
theory to account for arbitrary distributions through the positive-exponent velocity 
moments that are usually used to define the fiuid variables. This avoids the problem of 
divergent integrals. It also provides an obvious connection to the original fiuid criterion, 
and gives predictions that agree with experimental and numerical data. 
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